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Abstract 

We apply the method of principal value resummation of large momentum-dependent 

radiative corrections to the calculation of the Drell Yan cross section. We sum all next- 

0\ • to-leading logarithms and provide numerical results for the resummed exponent and the 
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corresponding hard scattering function. 



1 Introduction 

All-order resummation techniques have been in existence for some time now |]l], |^, mainly 
in an attempt to control perturbative QCD corrections to specific hadron-hadron scattering 
cross sections which are numerically large at a fixed order of perturbation theory p, Q. 
Such resummation formulas P| typically organize the large perturbative corrections in an 
exponential containing integrals of the QCD running coupling ^^(/i^) over various momentum 
scales and suffer from divergences associated with the Landau pole of as(/i^). 

In it was shown, in the context of the inclusive dilepton cross section, how to define 
a consistent resummation formula by a principal value prescription for the infrared (IR) 
coupling constant singularities of the exponent. This defines unambiguously the resummed 
perturbative series, as the exact exponent may be approximated by a finite series in the fixed 
coupling constant, thus avoiding an explicit and numerically important presence of higher- 
twist contributions through an arbitrary IR cut-off . The principal value prescription then 
defines a best perturbative resummation in the following sense: A perturbative exponent can 
be constructed as an asymptotic approximation to the principal value exponent and its degree 
(in as{Q'^)) is precisely determinable numerically as a function of Q^ as well. Thus, the often 
quoted "asymptotic nature of the QCD perturbative series" is given a precise meaning. 

In Q it was also shown that all large momentum- dependent perturbative corrections 
were contained in leading and next-to-leading exponents {El and Ej^l respectively). E^ was 
studied in that work, both analytically and numerically. -Ea^l is, however, also important for 
a precise perturbative description of the process and a full study of this exponent is necessary 
before the non-perturbative parton fiux is convoluted with the resummed hard part in order 
to yield the theoretical prediction for the cross section. 

The purpose of this work is to address these issues and perform a thorough analytical and 
numerical study of the exponentiated hard part, including all large momentum-dependent 
perturbative corrections. In a subsequent paper, we will use this hard part to calculate the 
dilepton cross sections at both fixed-target and collider energies and compare these with 
experimental results. 

This paper is organized as follows: in section 2 we briefiy summarize the new resummation 
formula and the expressions for El, both as an exact function of Q^ and as an asymptotic 
perturbative series, i.e., as a polynomial in as{Q'^)- This makes contact with earlier work 
1^ and defines notation and E^l as well. In section 3 we calculate E^l, in the spirit of 
0, both exactly and as an asymptotic perturbative series. The results are presented in 
both moment and momentum space, using a closed formula for the inverse Mellin transform, 
which has been derived elsewhere [Q . In section 4 we study the higher- twist inherent in the 
exact exponent of the resummation formula. We establish, in moment space, the asymptotic 
properties of the exponent in the non-perturbative region. These properties make evident 
the smooth behavior of the cross section for these extremely large moments. In section 5 we 
present numerical results for El and E^l in the perturbative and non-perturbative region, 
along with the corresponding perturbative approximations. We also present the hard part in 
momentum space, for a variety of values of Q^, including all large perturbative corrections. 
Finally in section 6 we summarize our findings and anticipate future, phenomenologically 
oriented work. Some technical details are given in an appendix. 



2 Principal- Value Resummation 

In this section we briefly summarize the principal-value resummation formula as applied to 
the dilepton production hard scattering function. Let us define kinematic variables for this 
reaction, denoted by 

hiipi) + h2{p2) ^ mn + X . (1) 

We denote s = (pi +^2)^? t = Q^/s, z = r/xaXb where Xa, Xb are the momentum fractions of 
the partons participating in the hard scattering. We shall concentrate on the diagonal-flavor 
quark- ant iquark hard parts, which alone are singular in the soft-gluon limit z — ^ 1, at a fixed 
order of perturbation theory. 

For quark-antiquark scattering, we obtain the principal-value resummation formula in 
moment space 0, §|, §[ 

Cjqq{n, a) = A{a)I{n, a), (2) 

where the function A{a) will be fully analyzed in a future paper, and where a = as{Q'^)/'n:- 
All n-dependence in ojqq is contained in the exponential 

/(n, a) = exp[E{n, a)L + E{n, a)NL] ■ (3) 

In eq. (Q), E{n,a)L includes the one-loop running coupling effects, and E{n,a)NL the two- 
loop effects. In p, ^ it was shown that, if we define the range of validity of perturbation 
theory by as(Q^)lnn < 1, this formula includes all large perturbative corrections. In 
particular. El resums all corrections containing more powers of Inn than of a, while Enl 
completes the resummation of corrections having equal powers of these two quantities. As 
usual, conclusions in momentum space may be drawn using the correspondence [|, ^ ra <-^ 
1/(1 — z). For El in the Deep-Inelastic Scattering (DIS) scheme, we can write the expressions 
i 

(4) 

(5) 
(6) 



E{n,a)L = a{g['h^-gi'h,), 




with 




h = h{n,t)=2I{n,t/2)-I{n,t), 




where 




/(n,t) =t|^c;c(^"^^~^) ln(l + (1/t) ln(l - 
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and where we define 

a02 \A^/ 





(7) 



Here g[ , g2 are numerical coefficients [^], while 62 is the first coefficient of the QCD beta- 
function. The explicit expressions for these coefficients are given in section 5. -Eatl is given 
by a similar, albeit more complicated, expression, which we will also give below. 



^A precise definition of tliis range will be discussed below, and in section 4. 



The principal value prescription makes the integrals in eqs. (@)-(0) well-defined at 
— jln(l — (^) = 1. P stands for an averaged mirror-symmetric contour connecting and 
1 above and below the real axis. It will be useful to consider two such special contours 
for analytical and/or numerical results. The definition of P, along with these two specific 
choices, appear in fig. 1. 

In P], a series expression for the integrals in eqs. (1)- (|^) was found, using the principal 
value contour proper, V: 



m=l 



TJllTTl 



, (1) -^ (^ ~^)mc^ 4-\ 

_^^(i) g 0:^_^s^^t), (9) 

where (a)^ = r(a + ^/^{(^), T being the Euler Gamma function, with 

£(x) = xe-^'Eiix), (10) 

and where the exponential integral is defined by the principal value 

Ei{x) =r r dyey/y. (11) 

J ~oo 

A perturbative version of eq. (|^) may be obtained by performing an asymptotic expansion 
of the special functions S{mt),S{mt/2) for mt,mt/2 ^ 1 in terms of In-' n/f^ and in a 
range where the Inn behavior is suppressed for all terms in the sums[§. For example, the 
second sum in (P) can be written as 






' ^ 777 '777 T^ ^r-rtlmP^'^ 



where use had been made of eq. (3.2) of 



Applying eq. (3.8) of 






E ^V^ = (-1)' E Cp+l-,^^ In^ n , (14) 



to eq. (0), one then obtains 



«^!'^ E %n;;^^M =«^f'^E^(-ir'E^P+2-,^^ln^n^ (15) 



From eq. (]T5|), we see that for n ^ 1, the suppression of Inn by 1/t can be quantified for 
a given power l/f^ by constraining the corresponding monomial In-* n, j = p. This leads 
to 062 Inn < 1 or, equivalently, ne~* < 1, as the range of validity of perturbation theory. 
We will further discuss this argument from a different point of view, when we investigate 
the non-perturbative regime in section 4. Of course, in eq. (|^) the sums have two different 
scales, the smaller one being t/2. Thus, if we want to define the perturbative regime for 
all three sums, we should restrict ourselves to n for which 2ab2lT^n < 1. In summary, the 
perturbative range in n is 

t 1 Q 

l<lnn<- = — — .^l<n<-f. (16) 

2 2ab2 A 

For n in this range, a large-t asymptotic approximation leads to the perturbative formula 

2 Af/'+l p+1 Afl'+l p 

E{n, a)L ^ E{n, a, iV)L = E E «' E 4p(^ ^^' ^ + E «' E 4(3) 1^^' ^^ (17) 

i=l p=l j=0 p=l i=0 

where the N^ determine the best asymptotic approximation, and where the coefficients are 
given by 

4(1) = -g^M-\-ir^^-^^2^Cp^.-,. 

4(3) = 9i'^b^2-\~ir'^-^^cp^r (18) 



The coefficients Ck are given by r(l + z) = YlT=Q^kZ^- Notice that El in eq. (|T^ becomes 
a polynomial in Inn whose precise powers are determined by the above large-t asymptotic 
approximation (perturbative approximation). Hence, the optimum number of resummed 
perturbative corrections, N = {A^/"}, is also determined in principal- value resummation. 
Numerical results for those numbers were given in [^ and will be incorporated and enlarged 
upon in sec. 5 below. 

3 Evaluation of the Next-to-Leading exponent in the 
perturbative regime 

Following 0, we may write the non- leading exponent E^l from eq. (0) as an integral over 
the function 



NL( rx^2iN (1) 2/, /, ^ ln(l + a&2 In A) (2) 2 

g, {a[XQ]) = -g, "^ ^^''^' \l + a^lnXf ^ '^ "" '{l + ab.lnXy 



+ ^ra'7^^-r^T?' ^ = 1'2. (19) 
l^i -t- au2 m Aj" „ „ 

In these terms, E^l is 



E{n, a)r,L = -l dc{^^^Y^) {/ Y^d^'i^ii^ " 0(1 " y)Q']) + 9^'Hi^ - OQ'])] ■ 

(20) 



From eqs. (|T^), (^) we may write the next-to-leading exponent as 

E{n, a)NL = a{g?Ji - g^^^ J2) + a\gf^ K, - g^K^), (21) 

with 

^--//<n^) (i + (iA)in(i-c)r ('') 

These integrals may be readily computed by separating the moment variable dependence 
from the coupling constant dependence, as in 0, and using the contour V of fig. 1. Care 
must be taken, however, regarding the singularity structure of the integrals over the various 
pieces of the contour. The semicircle contributions will now contain singular pieces of the 
form W6/6^, i,j = 0,1,.. which will cancel with similar singularities coming from the 
straight-line contributions. Let us consider, for example, the integral J2. Using the binomial 
expansion, we write 

J2 = E ^^^jr, (26) 



m=l ^'■ 



with 






,_,ln(l + (l/t)lnC ) 
(l/t)lnC)^ 

The semicircle contributions each contain an imaginary part, which cancels upon addition, 
and a singular real part: 



with (t = e~*. The straight-line integrations, on the other hand, give a combined contribu- 
tion: 



2 ' ' ' 2 ' ' . V o, ,, ^t ^ ^^^^ ■ ^^^^ 
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Upon addition we obtain the finite result 



— / dxe''\n^{\x\/{mt))+mV / dx—] . (29) 



= {h/hl){l-£{mt))-{h/{2hl))K{mt), (30) 



with 



A(x) = xe-^(7r2 - / dyenn^ {\y\ / x)) . 



(31) 



The function A(x) is expressible in terms of Incomplete Gamma functions, as we will show 
later in this section. 

Let us first give the results for the integrals Ji, Ki, K2 appearing in eqs. (p2D-(p5D. 
Performing the ^-integration for the integral Ji and expanding in terms of the "Pochhammer 
symbol" (1 — n)m as before, we have 



^ &3 f^ (1 - n)m 
C'2 m=l "^' 



where we define 



and 



Hence we obtain 



jp 1 



J- (t/2) - J™ (t) + J- (V2) - JU^) 
iln(l + (l/t)lnC) 



+ (l/t)lnC 



2^''^^'^^ 



JU^) - /p^cr-^YTTiW^ = ^^^"^'^ • 



^^_&3 y. (1-nU jl 



L/t)ln(^ m 



A(mt/2) - A(mt) I + £{mt/2) - £(mt) 



(32) 

(33) 
(34) 

(35) 



Working similarly for the integrals Ki, K2, eqs. (5.4), (5.5), we find 
and 



ah 



1 £(L_^{^(^i/2)-^(^t)}, 



'9 1 

^ m=l 



m\m 



1 00 d _ r,^ 

^2 = 4 E ^^^{^("^t) - 1} . 

Putting everything together we find for the next-to-leading exponent the expression 



(36) 



(37) 



Ein, a 



NL 



\ m=l 

1-n), 



(1 -n), 



m\m 



S{mt/2) - S{mt) + -(A(mt/2) - A(mt)) 



+9^2'^ E 



m=l 



Till 



S{mt) -1 + -A(mt) 



+W&2) -^r^E 

I m=l 

(2) v^ (1 — "')m 



(2) v^ (1 ~ fT'Jm 



ni\ni 



£{mt/2) - £{mt) 



-9^'^ E 



m=l 



ml 



£{mt) - 1 



(3^ 



To recover the resummed perturbative series, we have to perform asymptotic expansions on 
the functions £{x), and A(x), eq. (^Tj). 



As with the function S{x) = xe ^Ei{x) (see eq. (|ri|)), we will now derive an asymptotic 
formula for the function A(x), eq. (|3ll) . The definition of this function is 

A(x) = xe"^(7r^ - / (i?/e^ln^(|y|/x)) ~ -xe"^|72(x) - 21na;7i(a;) + In^ x7o(x) i , (39) 

where 

^n{x)= / dyeHn^y, (40) 

Jo 

the difference between the two expressions in eq. ( ^91) being of order (In £)^- The functions 
7„(x) may be expressed in terms of Incomplete Gamma functions |T0| 

7„(x) = J_im (^-j J^ dye^y' = £m (^-j x'+^T{t + l)7*(e + 1, -x) . (41) 

Using the asymptotic expansion for the Incomplete Gamma function, or simply integrating 
eq. (|40|) by parts, we obtain 

-"•^-^-(DWJl+lti).}, ,42) 

where N' will be defined shortly. Realizing that (— e)p is a polynomial of p degree in e, with 
roots ej = j, j = 0, ...p — 1, we finally obtain 



7o(a;) ~ e"'. 



liix) 
72(3;) 



1 ^'"1 p! 
In a; T^ — 

X p=0 -^ 

1^2 ^ _ 21nx ^^1 p! ^ 2 ^^1 [vl/(p + 1) + 7]p! 



p=o •^' -^ p=l ^ 



(43) 



Putting everything together we obtain the asymptotic expression 



AW._2i:Iik±Jl±^, ,44) 

p=l "'' 

where \E' is the logarithmic derivative of the F-function. As with A^ = {^i^} iii ^q. {^), 
N' depends on x and different A^"s will be determined for the different sums in eq. (^). 
A similar method was introduced in |^ for E^. Let us separate, therefore, the sums in 
eq. (^) and derive optimum numbers of perturbative terms by direct optimization of each 
sum separately. Using eq. (^) along with the corresponding asymptotic formula for S{x), 
0, we obtain the perturbative asymptotic approximation: 

jyJVL 

AjNL 
2 I 00 / -1 \ 



p=0 



tP j~^ m\mP~^^ 



N^^ 



i<yh/bl)g['^T.^-^mP+l)+l]t^^ "^ 



p=l m=l 



,2^ Ji)^ p!r.T,.„ , .^ , .1^ (l-^)r 






p=0 t' ^1 ^^-^^ 



m=l '"'■ 



-M3/&^k^^^ES[^(p + i) + 7]E^^^ 



m 



p=l ^'^ ^1 "^'"^'^ 



(„/y,p."£4?!£(-" 



p=o ^^ ^1 "^!"^^+' 
(2) ^' p! ^ (1 - n). 



-(«/&2k2 E^E^^^y;;^ 

p=0 rn,=l 

+ («/&2)^f E ^M^ (45) 

The set of numbers A^' = {N^^}, i = 1, ...6, is determined by direct optimization of 
eqs. (pSD, (^51), as in the leading case. Then, performing a Stirhng approximation on the 
rz-dependent summations we observe that the fifth and sixth sums contain at most one less 
power of In n than of coupling constant and hence we drop these for consistency. Therefore 
we end up with the perturbative formula 

E{n,a,N')^L = Y: E a^t^sf^l^it)^ n, (46) 

i=l p=2 j=0 

where the perturbative coefficients are 



J! 



sSi2) = [g['\bs/bi)-g['^/hn-\-ir^ 



,(P-1)! 



Cn 



r. ^'~'' 



sf,^{S) = 9['\b3/bln-\-ir'^-^j^2'-'[^{p)+lK-, 
C(4) = -9?\b,/blM-\-ir^^-P^[^{p)+j]c,^,. 



-- - Jl 



(47) 



The perturbative exponent in moment space, resumming all large perturbative corrections, 
becomes, after eqs. (0), (|T|), (||), (0) : 

E{n, a, ^f) = E{n, a, N)l + E{n, a, N')nl, (48) 

with Af = Af{t) = {N;N'} = {N[, N^, N^; N^^^, N^^, N^^, N^^}, all numbers being 
functions of t obeying the relations 

N^{t) = N^{t/2), iVf ^(t) = <^(t/2) = N,\t), iV3^^(t) = iVf ^(t/2) , (49) 

which give a total of three independent functions. The resummed hard parts in momentum 
space can now be obtained, given the exponent eq. (|48|) . In the perturbative regime, a closed 
formula for the hard parts giving all large perturbative corrections is [§: 



^(t^.^'AT) ^ 



7T 1-z) 



r (l + Pi (737.".-^)) sin (vrPi (y3^'«'-^)) 



(50) 



Pi (n, a, AT) = ^^E{n, a, AT) . (51) 

omn 



I{z,a) = 5{l-z) 
with 



The range of validity of eq. (|50D is 2ah2 ln(l/(l — z)) <1 and in fact, numerical calculations 
in P] suggest that this upper limit is in good agreement with the peak position of the exact 
leading exponent, 062 ln(l/(l — zl^)) ~ 0.6. For larger 1/(1 — z), the higher-twist inherent in 
the resummation formula, takes over, causing the perturbative approximation to break down. 
In the next section we will determine the behavior of the higher twist analytically, as a large- 
n asymptotic approximation in a regime complementary to the perturbative one, namely for 
a62ln(l/(l — z)) > 1. In section 5 we will present numerical results for E{n, 0)1/^1 and 
E{n, a, N/N')l/nl along with the sets of numbers N and A^', for various values of a. 

4 Behavior of the exponent in the non-perturbative 
regime 

Before proceeding to numerical results, it is important to understand analytically the behav- 
ior of the exponent for asymptotically large values of n -beyond the perturbative regime. This 
will reveal both the behavior of the higher twist inherent in our resummation formalism, as 
well as that of the cross section itself at the very edge of phase space z ^ 1, or equivalently, 
n — i> 00. It will also confirm and extend the corresponding numerical calculations, which are 
quite sensitive at this range. In |^ a heuristic argument for the behavior of E{n,a)L was 
given, which was verified by numerical results showing that 

lim E(n,a)L = —00 . 

Here we will elaborate on this argument, developing asymptotic formulas for both E{n, a)^ 
and E{n,a)jsfL in the non-perturbative range 1 < 062 Inn < cx3. These formulas will be 
checked against numerical results in section 5. 

10 



Let us begin with the leading exponent, eq. (|^). We will concentrate on the integral 
I{n,t), eq. @. After a change of variables C ^ 1 — C; it reads 

Iin,t)=t f ^iy(n,C)ln(l + (l/t)lnC), (52) 

with the n-dependent weight given by 

W{n, C) = [(1 - C)"^' - 1] = e"^""'^ '"(^) - 1 . (53) 

The branch-point singularity of the above integrand is at Ct = e~*- On the other hand, the 
weight approaches the value —1 when| 

( > 1 — e"^^^ = no(n) . (54) 

As n varies, we may have either (t < '^o(^) or (t > no{n). The former inequality is satisfied 
whenever 

n<l + —^^ = Ciit). (55) 



In fact, for reasonably small values of a, Ci(^) — Ct"^ = e* = Q^/A^ and hence eq. (|55| ) is 
equivalent to 062 Inn < 1, which we recognize as our perturbative regime for this particular 
integral. As before, some of the integrals in El and E^l contain the scale t/2, but again 
we take the intersection of the regions defined by Ci(^)) Ci(V2) to define the narrower non- 
perturbative regime n > Ci{t) ^ e*. This leaves a rather uninteresting intermediate region 
Q/A < n < Q^/A^, where various scales impose different asymptotic behaviors to different 
pieces of the exponent. 

Therefore the non-perturbative asymptotics are described by the conditions 

n > Ci(t) ^ e* = ^, W{n, () ~ -1, C > no{n) ^ -, (56) 

where the asymptotic behavior of /(n, t) is given by 

I{n,t)^c^-tf ^ln(l + (l/t)lnC), (57) 

with Vo the principal- value contour proper, fig. lb, but starting fron no(n), rather than 0. 
We have computed the integral on this contour in the appendix. The result is 

/(ra,t)> c^tl -tL{n,t) In L{n,t) + ln(l/rao(n))| =tG{n,t), (58) 

n> (lit) = 1+ . \ ^ -e* . (59) 

Inf 



^This value is the source of divergent behavior as rt ^ cx). 
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Let us now digress for a while, to complete the discussion in section 2 regarding the 
definition of the perturbative regime. We will confine our discussion to the integral I{n,t) 
as given by eq. (|^). In section 2, we showed that a general expression for this integral is 

I(n, t) = - V il^4^^(mt) . (60) 

We also showed that for 

t > 1, K n < e*, (61) 

we obtain the perturbative approximation 

Hn, t)c^-Y: Tp^-^r' E ^P^^-j—r 1^' ^ = ^(^' ^' ^) • (62) 

p=0 '' j=o J- 



Even though the second condition in eq. ( |6lD is the one which suppresses the logarithms 
and hence makes an asymptotic expansion possible in I{n,t,N), only the first condition in 
eq. (|6T|) was used to derive this approximation. Hence, the question might arise whether 
the perturbative approximation (^) is valid independently of the second condition in (pl|), 
which defines the perturbative regime in n. We now see from eqs. (|^), (0), (0), however, 
that for n > Ci{t) ~ e*, the integral I{n, t) is approximated by I{n, t)> which is clearly a non- 
perturbative expression. Furthermore, as we will see in the rest of this section, expressions 
like I{n,t,N) and I{n,t)y are numerically very different from one another, for n > e*. | 
Thus, when t > 1, the perturbative expression I{n,t,N) approximates I{n,t) for n < e* 
while the non-perturbative expression I{n,t)y approximates I{n,t) for n > e*. This justifies 
the characterization of the corresponding n-ranges as perturbative and non-perturbative, 
respectively. 

Returning now to our calculation, we find for l2{n,t), eq. (0), 

l2in,t)y^t\nL{n,t)=H{n,t), (63) 

with L{n,t) given by eq. ([591) . Our leading exponent is then given, in this region, by 

E{n,a)L> = ag['\2I{n,t/2)^-I{n,t)^]-agi'h2{n,t)^ . (64) 

From eqs. ( ^^ - |6^) we find that, when n > e*^*, 

'lnn\ 



and hence 



G{n,t) ~ —Inn In 



(1) (1) /] \ 

E(n,a)L>~-^ln21nn-^ln — . (65) 

02 02 \ t I 



The first term above dominates the behavior as a logarithmic divergence in n, with a negative 
coefficient. This point was also made in P, but the precise asymptotic behavior, eqs. (|58D - 
(|63D, valid in almost all of the non-perturbative range n > Ci{t), was not given. 

^I{n,t)y, on the other hand, is not even defined for n < e*. 
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We can similarly find the large-n asymptotics of the next-to- leading exponent, 

E{n,a)ML = a[g?J^{n,t)-g'i\h{n,t)]+a^[gf^K^{n,t)-gf^K2{n,t)] . (66) 

Concentrating on the integrals Ji and i^i, eq. (pi]), and using the results of the appendix, 
we can write 



Ji(n, t)> ~ ^{i^(n, t/2) - F(n, t) + H{n, t/2) - H{n, t)|, (67) 



b. 

with H{n,t) defined in eq. (|63D , where 

t 



F{n,t)=^-{\n'{L{n,t))-n'y (68) 



and 



Ki(n, t)> ^ -tlH{n, t/2) - H{n, t)\ . (69) 

Similarly, for J2 and K2 we find 



J2(n,t)>^-^P(n,t), (70) 

"2 



where 



and 

i^2(n,t)> ~ ^{ w^ + 1} = Qin,t) . (72) 

Notice from the above expressions that the quantities J2, K2, which were subleading in the 
perturbative regime, remain subleading in the non-perturbative regime as well. In fact, the 
dominant large-n behavior in E{n, 0)^1 comes from the function F{n, t), eq. (pSI). For values 



n > 1 + — — \ y ~ e*[i+<^^l , (73) 

the next-to-leading exponent is dominated by 

E(n,a).,>^a,i^^|-ln (^-2!- [—) \ ^ '^^'4^,'^ [—) ' (^4) 

Therefore, E{n,a)NL> also tends to —00, though more slowly than E{n,a)L^. This only 
happens, however, for huge values of n, eq. (|73D , quite inaccessible even to numerical cal- 
culations. As we shall see in section 5, within the numerically accessible range of n, the 
singularity manifests itself through the presence of vr^, eq. (|68D , which dominates, making 
E[n, a)NL look like it is tending to -|-oo within that range. This is a demonstration of the 
importance of the above asymptotic formulas where numerical calculations break down. 
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Numerical results for the Exponent and the Hard 
Part 



In this section we present numerical results for the exponent of the resummation formula 
as well as the qq hard part, at various values of the kinematic variables, and including the 
large perturbative corrections. In our opinion, these results exhaust our knowledge of the 
perturbative phase of the theory. The various constants appearing throughout this work, are 



9i 

&2 



(1) 



;iica 



9'^' 



3^ 



(2) 
^1 



Cf 



c, 



67 
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9 



2nj)/12, 63 = (34C2 - (lOC^ + QCF)nf)/AS. 



(75) 



For the rest of this paper, we will use the values ||^ A = 0.2GeV, Uf 
5, 10, 90GeV, which correspond to a = 0.075, 0.061, 0.039 respectively. 



4 and Q 



5.1 The resummed exponent 

Let us first present numerical results for the exponent in moment space. Approximate 
conclusions in momentum space may be reached through the correspondence n ^^ 1/{1 — z). 
As in [^, let us first give the optimum numbers of asymptotic terms for eqs. (^5D, (^), 
for the above values of Q, computed at the intermediate value n = 30.0 As we have noted 
before, the set A/" depends on t = 1/(062) = ln((5^/A^). 

TABLE 1 



Optimum numbers of asymptotic terms for E{n, a)L 


E{n,a)NL 


ait) Nt N^ N^ ATf^ N^"^^ N^"^ 


N^' 


0.075(6.4) 15 5 1 5 1 


5 


0.061(7.8) 2 6 6 2 6 1 


6 


0.039(12.2) 4 12 11 4 11 4 


11 



Given these values we may compute the perturbative exponent as given in eqs. (|48|) , {^), 
(^), ( P7D and compare it with the exact exponent E{n,a) = E{n,a)L + E{n,a)]\fL. 
The latter is computed numerically on the box contour [^ P, fig. Ic. In fig. 2 we have 
plotted the result of these comparisons for the above values of a. For completeness we have 
also included comparisons of the leading exponents alone. Notice that we expect the peak of 
the exact exponents to grow as a gets smaller. This is due to the fact that the higher-twist 
tends to suppress the exponent, hence, the smaller the a, the less the suppression. Notice 
also that the exact and perturbative curves are in excellent agreement near the peak of the 
exact curves. We will comment on this in some detail later. Once past the peak of E{n, a), 
the higher- twist completely dominates the exponent. 

Using the results of sec. 4 we may compare the exact exponents with their large-n 
asymptotic expressions, eqs. (^), (ff^ in the non-perturbative regime. This is done in fig. 
3 for various values of a, within the maximum range of n attainable before we run into 



"* As with the leading exponent, there is very httle dependence of these numbers on n, when the latter 
is well within the perturbative regime. Strictly speaking, for Q = 5GeV, n = 30 reaches the limit Q/A for 



N{^, Ni 



NL 



N. 



NL 



However, these numbers are unchanged for lower values of n, as well. 
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round-off errors for the exact exponents. Note that the agreement between the exact curves 
and the analytical approximations improves as n increases. 

In fig. 4a we have plotted separately E{n, a)NL and the corresponding perturbative 
approximation for various values of Q. We see that both curves have the same shape for 
all values of Q used, but that at low Q-values numerical agreement is only fair. We can 
trace this to the determination of A^' and A^^^ in particular, see eqs. (pSf), (0), and table 1. 
Using the exact expression for A(a;) as given by the first equality of eq. (|39D , we see that the 
TT^-factor is obviously of higher-twist origin and can not be expressed as a polynomial in a. 
However, its effect on the third sum of eq. (^) may be isolated. The equation determining 
N^^ may be written 

The LHS gives, for Q = lOGeV, the values {1.798,0.363} for n = {10,30} respectively. At 
the same time the best approximation on the RHS is given for A"^^ = 1 and equals the 
minimum possible values {1.219,2.212} respectively, which is admittedly not a very good 
perturbative approximation, with an error relative to the perturbative value of {47%, 83%} 
respectively. If, on the other hand, we subtract the effect of the vr^-term in the LHS of 
eq. (^) that enters the definition of A(x) in eq. (0), which may be computed analytically 
to be 

^[(l-e-*/^r^-l], 

we obtain the values {5.008,8.923}. These can be approximated by the perturbative ex- 
pression for A^^^ = 3 — 4 yielding the values {5.568,8.395}, in very good agreement with 
the previous ones, with an error of only {10%, 6%}. For values of Q near the Z-mass, on 
the other hand, numerical agreement between the two curves is excellent in the perturbative 
regime. To summarize, the accuracy of the perturbative approximation for E{n, ojnl is fair 
for low Q because the higher-twist is much more important, and is suppressed more slowly 
with increasing Q, than in the integrals giving the leading exponent, E{n, a)^. 

On the other hand, in fig. 4b we have plotted separately E{n, a)NL in the nonperturbative 
regime and have compared it with the corresponding asymptotic approximation E{n, a)^^^, 
up to the maximum values of n where a reliable numerical calculation of the former is 
possible. Notice the excellent agreement as well as the fact that these quantities have a 
positive increasing behavior in this range, something exhibited by the dominance of the 
TT^-term in eq. (p8|), as remarked in sec. 4. Of course, for even larger values of n, the 
next-to-leading exponent tends to — oo as suggested by eqs. (0), (|7^). 

Let us now return to fig. 2 and discuss in more detail the peak positions of the exact 
exponents. In it was found that the peak of E{n, a)^ is pretty close to the perturbative 
boundary n ^ Q/A. Adding E{n, a)NL decreases both the magnitude and the position of 
the peak, as seen in the above figures. Let us therefore define some characteristic points 
in these curves which will describe, within a controlled arbitrariness, the interface between 



^ These values are inaccesible to an exact numerical calculation but are very easily handled numerically 
through the approximate expressions, eqs. (n6[|72). 
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the perturbative and non-perturbative regime and a corresponding theoretical uncertainty 
in the resummed perturbative cross section. One such point that we have aheady mentioned 
several times, e*'^ = Q/A, separates the perturbative from the exact analytical expression 
for the exponent. This, from an analytical point of view, is clearly the boundary between 
the perturbative and non-perturbative regimes. A second point, already mentioned in 0, is 
the peak position of the exponent, Cp '■ E{(p, a) = max. In fact we may define separately 
the peak positions for the leading exponent]^, (^ and the full one, Cp- The table that follows 
gives these numbers as a function of Q. 





T/ 


^BLE 2 




Boundaries in n - 


- space parametrizing 


the perturbative/nonperturbative interface 


QiGeV) 


g/A 




Cp 


5 


25 


37 


14 


10 


50 


100 


37 


90 


450 


1875 


570 



Notice in the above table that inclusion of E^l creates good agreement between the order 
of Q/A and CpiQ)- That is, while 2ab2 In (^ is close to unity with good accuracy, ranging from 
1.12 [Q = 5GeV) to 1.24 {Q = 90GeV), inclusion of the next-to-leading exponent moves the 
peak so that 20:62 In Cp(5; 10; 90GeV) = 0.82; 0.92; 1.03, is close to 1. In fact these results 
strongly support the view that the peak of the exponent occurs at the limit of validity of 
perturbation theory and approximately coincides with the analytical result, namely Q/A. 
The fact that we have two different scales in the problem, t and t/2, is reflected in the minor 
differences between Q/A and (p in table 2. This difference may be used to provide a small 
theoretical uncertainty in the resummed perturbative hard part but, for the purpose of this 
work we will consider e*/^ = Q/A to be the boundary of the perturbative regime in moment 
space. In momentum space, this converts into a perturbative regime given by 

< z < 1 - A , (77) 

and this is what we will use for deflning all resummed perturbative quantities below. 

5.2 The resummed hard part 

A formula for the hard part of the qq cross section, which certainly contains all the pertur- 
bative corrections that eq. (^) does, but where the higher-twist is not separated, may be 
obtained by "extending"[] eq. (^Dp throughout the whole range of z, by numerically calculat- 
ing both the full exponent E{1/{1 — z),a) and its derivative 



z' J d\n{l/{l-z)) \l-z 



^This point was denoted by rii in [pj. 

^ Strictly speaking, the inversion of the MeUin transform [pj, was performed only for the perturbative 
exponent. The corresponding formula for the hard parts, involving power-counting arguments for all large 
perturbative corrections, simplifies into eq. (pOl) only for the polynomial approximation of the exponent as 
powers of a and Inn. 
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on the "box" contour P, fig. 1: 



I{z) = 6{1 - z) 



-.Ei^ 



y^r(i + A(^,„)).i„(.p,(^.„ 



(7f 



As is discussed in [Q, in the physical cross section we may perform an integration by parts 
and explicitly get rid of the plus-distributions as well as the (5-function in the above formula. 
Then, the relevant momentum-dependent quantity contributing directly to the cross section 
as the perturbative hard part is 



E 



On the other hand, from the previous section, and using the results in tables 1 and 2, we 
can now find the resummed perturbative hard part by separating the higher twist. Since the 
latter has the effect of reducing the exponent till, as n — ^ oo, it goes to — oo, we may remove 
its effect by setting the perturbative hard part equal to zero in the non-perturbative regime 
1 — A/Q < z < 1. We thus arrive at the following perturbative formula: 

n{z, aM)= (i - ^ - ^) //^'4(rz^^ (^ + ^^ (t^' "'■^)) 

xsin(7rPi(^,a,Ar)). (80) 



To avoid repetition, let us define the notation that, for any resummed quantity of interest 
i?, R{N)l stands for this quantity computed using eq. (|80D with E{N)l only, R{M) using 
the same equation but with the full E{Af)) and similarly for Rl and R, but using eq. ( [7D[ ) 
instead. We will now present numerically the hard part computed according to these various 
definitions. 

In fig. 5 we have plotted the functions H{z,a)L, H{z,a) as well as their perturbative 
counterparts. Some comments on these curves are in order. As a general observation, 
the hard part calculated with the total exponent agrees very well with its perturbative 
approximation in almost all of the z-range and the agreement improves with increasing Q. 
In particular, near the Z-mass the agreement is excellent. In the neighbourhood of the 
non-perturbative regime, z > 1 — A/Q, we observe some interesting behavior, naturally 
more pronounced at low Q-values. Let us first concentrate on TC{z,a,J\f)L, computed from 
eq. (^Of). Going for simplicity to n-space, we observe that for any value of Q, the polynomial 
Pi{n,a,Af) is a positive increasing function of n in the perturbative regime but may be 
of order unity for sufficiently large values of Q. This was actually anticipated in P], and 
does not contradict the power-counting of a and Inn developed in that reference, which 
established that Pi contains at most equal powers of these quantities and hence is anticipated 
to be < 1 in the range 2a62liin < 1. In fact this expectation is numerically accurate 
for the normalized series Pi/Qi , for all values of Q. The somewhat large (and arbitrary 
from a resummation point of view) value of the color factor g^ = 2Cf = 8/3, which is 
automatically included in the resummation as an overall multiplicative factor of E{n,a, Af)L, 
as it should, can make Pi somewhat larger. To be concrete, for Q = 5, 10, 90GeV, Pi(n = 
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Q/A)l = 0.79, 1.67, 2.74 respectively while the corresponding normalized series has the 
values 0.30, 0.62, 1.02, in agreement with the power-counting estimate. For the values 
Pi(n, a,A/')L > 1 at corresponding values of Q, the sinusoidal dependence of the hard part 
will produce a decreasing behavior which may even become oscillatory for sufficiently large 
Q. The decreasing behavior is shown in fig. 5b, whereas the oscillatory behavior occurs for 
Q = 90GeV but is not shown. These effects are obviously concentrated in the neighborhood 
of 2; ~ 1 — A/Q in momentum space, but are still part of the perturbative regime due 
primarily to the overall color factor. 



Similar remarks apply to the "exact" curves, produced through the formula ([79|), where 
the sign of the integrand will become negative whenever Pi > 1 or Pi < 0. We note that, at 
Q = 5GeV both P/" and Pi are less than 1 in the entire perturbative regime, at Q = lOGeV 
P^ > 1 for 15 < n < 38 while Pi < 1 in the entire perturbative regime and at Q = QOGeV 
P/"; Pi > 1 for 39 < n < 935 and 27 < n < 319 respectively. These ranges will create 
oscillations in the corresponding curves for the hard part, as is shown in fig. 5b for 7i{z, a)L- 
As we tend towards the peak position of the corresponding exponent, table 1, P/"; Pi will 
decrease again, change sign as they cross the peak position, and become negative in the non- 
perturbative regime. There is no danger from the poles of the F-function when Pi becomes a 
negative integer, because sin(7rPi), that multiplies it, cancels these poles. The corresponding 
curves will decrease thereafter, reaching a finite value at z = 1 since, as we have seen in sec. 
4, the exponent of the resummation formula ([79|) tends to —00 as z ^ 1. We indeed find 
numerically that, for Q = 5, 10, 90GeV, n{z ~ 1, a)L ^ -2.17; -10.85; -644.0 respectively, 
while T-C{z ~ l,a) ~ —0.72; —2.23; —151.4. | The corresponding curves for the cross section 
will be bounded for values of r near the edge of phase space r ~ 1. 

6 Conclusions 

In this paper we have completed the principal- value resummation of the Drell Yan cross 
section in the DIS scheme, by including the next-to-leading exponent to previously existing 
results 0. The resulting exponent E = E^ + E^^ includes the large perturbative cor- 
rections appearing as plus-distributions in the perturbative calculation of the cross section 
and produces a resummed hard part ^ that is finite throughout the entire kinematic range 
and independent of IR cutoffs. This principal-value exponent E may be approximated by 
a perturbative expression, P(A/'), obtained in moment space as an asymptotic series in the 
perturbative regime 1 << n < Q/A, and by a non-perturbative expression, E{n)-^ when 
n > Q^/A^. The former contains a finite number A/'(t) of perturbative terms, which depends 
on three independent functions of t = ln((5^/A^), while the latter tends to — 00 as n — > cx). 
This makes the hard part a finite function of z and the resummed cross section a bounded 
function of r. 

The accuracy of the perturbative approximation of the exponent improves with Q, as 
expected, being good at fixed-target values and excellent near the Z-mass. Ej^l{M) approx- 
imates fairly its exact version at low Q-values, due to large higher-twist contributions, and 



^The values quoted for Ti, are taken at the highest value of z, equal to 0.999999, before machine round-off 
errors are encountered by the program. 
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very well near the Z-mass. In short, E^l is numerically significant, as expected, in all its 
versions and all regions. 

The hard part for the Drell Yan cross section, 7i(2, a), whose general form was given in 
IP], is reproduced accurately in the perturbative regime by its perturbative approximation 
TC{z,a,Af) except near the edge z = 1 — A/Q, where both expressions peak and then 
oscillate (reaching negative values) due to higher-twist and/or large color factors. As before, 
agreement improves substantially with increasing Q, making the theoretical uncertainty in 
the resummation, if defined roughly as the difference between the above two quantities, a 
decreasing function of Q. 

We are now in a position to input this resummed hard part, with the properties de- 
scribed above, directly to a convolution with the non-perturbative parton flux for various 
experimental situations, to obtain the corresponding resummed cross section. Due to the 
multitude of experimental cases to be covered, the associated numerical complication of the 
corresponding calculations, the need for an accurate definition of the uncertainty in the re- 
summation and the essentially phenomenological character of these matters, we reserve them 
for a subsequent publication. 
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APPENDIX A Example of a large-n asymptotic calculation 

In this appendix, we will show in some detail how to obtain the asymptotic expressions 
of the various integrals in the exponent, in the non-perturbative regime, n > ^. As an 
example we will work out the integral Ji, eq. (^21) . Performing the y-integration, we have 



'1) 



Ji(n,t) = r^{ Ji;i{n,t/2) - Ji;i(n,t) + Ji;2(n,t/2) - Ji;2(n,t) 



where, after a change of variables C ^ 1 — C, 
and 



Jl;2{n,t) 



dC, 



+ (l/t)lnC 
1 



PC^^^'^^r+(i/t)inc- 



^2) 
^3) 



W{n, () is given by eq. (|53|). 

Focusing on the first of these functions we can write, after the discussion in sec. 4, its 
large-n approximation as 

rfCln(l + (l/t)lnC) 



Ji;i{n,t): 



4) 



ho C l + (l/t)lnC ' 

where the contour Vq extends from nQ{n) to 1 above and below the real axis and comprises 
a straight section and a semicircle around the singularity, with radius 6, see fig. la. 
The semicircle contributions to the integral are 



^</l;l(n) + ^Ji;i(U) 



_T3 r rO iSe'Ode ln(l+(l/f)lnCt + (<5ACt)e''') l 

= -Re{ttJ^de[\n{^)+te]} = -'4. 
On the other hand, the straight-line contributions are 



^5) 



-Ji;i(+ie) + -Ji;i(-ze) 



^^\Jno{n) X l+{l/t)lna; J kt+S x l+{l/t)lnx 



1-^4) + ^1-' 



t \no{n) 



In^ (4) . (86) 



From eqs. (]85|), (p6|), we obtain the desired asymptotic formula for the integral in the non- 
perturbative regime: 



Ji;i(n, a); 



In^ 



t \no{n) 



n'\^F{n,t). 



(87) 



Notice that the functional form of the divergent n-behavior in eq. (|87D , in terms of the 
quantity 

L{n,t)^hn(^]-1, 
t \no{n)J 

which contains the perturbative/nonperturbative interface, is dictated by the structure of 
the IR singularities in eq. (|86|) (which, of course, cancel). The rest of the integrals in sec. 4 
can be worked out in exactly the same way. 



20 



FIGURE CAPTIONS 

Figure 1. Contours for evaluating the Principal- Value exponent. 

(a) General definition. 

(b) Principal- Value contour proper, V. 

(c) "Box" contour P, for numerical evaluation. 

Figure 2. Exact and perturbative exponents. 

(a) Q = 5GeV: Dot=E{n,a)L', Dash=_E(n, q;,A/')l; 
Solid=_E(n, a); Dot-dash=_E(n, a, A/"). 

(b) Same as (a) but for Q = lOGeV. 

(c) Same as (a) but for Q = 90GeV. 

Figure 3. Exact and nonperturbative (large-n) exponents. 

(a) Q = 5GeV: Dot=E{n,a)L', T) ash.= E{n, a) l>', 
Solid=-E(n, a); Dot-dash=_E(n, «)>. 

(b) Same as (a) but for Q = lOGeV. 

(c) Same as (a) but for Q = 90GeV. 

Figure 4a. Exact and perturbative next-to-leading exponents. 

Solid=E(n, 0.075)7vl; Dot=E(n,0.075,A/')7VL; Short dash=E(n, 0.061)7vl; 
Long dash=E(n, 0.061, A/") iVL; Dot-short dash= E {n, .039) nl, 
Dot-long dash=E(n,0.039,A/')7VL. 

Figure 4b. Exact and non-perturbative (large-n) next-to-leading exponents. 
Notation as in 4a, but with the substitution E{n, a,J\f)j\fL -^ E{n, a)NL^. 

Figure 5. The resummed hard function 7i{z, a). 

(a) Q = 5GeV: Dot=7^(z, a)^; Dash=n{z,a,Af)L; 
So\id=T-C{z,a); Dot-dash=H{z,a,J\f). 

(b) Same as (a) but for Q = lOGeV. 

(c) Same as (a) but for Q = QOGeV. 
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